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Scalar electrodynamics can be used to investigate the formation of cosmic strings 
in the early universe. We present the results of lattice Monte Carlo simulations 
of an effective three-dimensional U(l)+Higgs theory that describes the equilib- 
rium properties of finite-temperature scalar electrodynamics near the transition. 
A gauge-invariant criterion for the existence of a vortex is used in measuring the 
properties of the vortex network in the equilibrium state both in the Coulomb and 
in the Higgs phase of the system. The naive definition of the vortex density be- 
comes meaningless in the continuum limit and special care is needed in extracting 
physical quantities. Numerical evidence for a physical discontinuity in the vortex 
density is given. 

The traditional picture for defect formation in gauge theories relies on the 
Kibble mechanismlll: the effective potential changes shape at the critical point 
and the field relaxes to different minima in space-time points separated by 
more than a correlation length. However, strictly speaking this approach can 
only be applied to transitions in which a symmetry breaks spontaneously, and 
gauge symmetries cannot be broken □. In some cases, e.g. in the electroweak 
theory, the transition line ends at a critical point after which there is only an 
analytical crossover between the phases a. It is clear that some changes must 
be made to the picture of defect formation in cases like this. 

We will here present the results of the simulations of three-dimensional 
scalar electrodynamics, i.e._the Abelian Higgs model. The details can be found 
in the original publications!!!! The theory describes the equilibrium properties 
of finite-temperature relativistic scalar electrodynamics near the transition 
and has the same form as the Ginzburg-Landau theory of superconductivity. 
We use lattice Monte Carlo simulations to study the behavior of thermally 
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generated vortices, i.e. line-like topological defects, in the two phases of the 
theory. Similar studies have previously been, carried out for the globally sym- 
metric theory using the Langevin equation u. Our final aim is to understand 
non-perturbatively the process of defect formation, but we find it necessary to 
understand the equilibrium behavior first. 

The continuum theory is given by the action 



S= d 



-;3 , 



(i) 



where Fij — diAj — djAi, Di = di + iAi and y is renormalized in the MS 
scheme at ji = 1. In this notation all quantities are dimcnsionless. If we write 
4>(x) = v(x) exp[i , y(x)], the action is invariant under the gauge transformation 

7(x)->[7(x)+0(x)], r> Ai(as) ^ Ai(x) - diB(x), (2) 

where [X] v = X + 2nn such that [X]^ G (— ir, ir}. The relations of x and y to 
the parameters of the 4D scalar electrodynamics and superconductors, as well 
as the determination of the phase diagram, have been discussed elsewhere u. 

The phase diagram of the system consists of the Coulomb and the Higgs 
phase. Although a mean-field analysis suggests a symmetry-breaking transi- 
tion, it turns out that there is, in fact, no local order parameter. The mass of 
the photon acts as a non-local order parameter, being non-zero in the Higgs 
phase and vanishing in the Coulomb phase At small x, the transition is of 
first order as predictecLby perturbation theory, but at some critical value of x 
it becomes continuous E3 (See Fig. [l^,). 

For numerical simulations, the theory must be defined on a lattice. We 
use the non-compact formulation, which means that there is a real num- 
ber cti(x) corresponding to the continuum gauge field Ai(x) on each link 
(x, x + i) between the lattice sites. On each site there is a scalar field 4>{x) = 
v(x) exp[ij(x)]. The lattice analogue of the gauge transformation (||) is 

j(x) ->• [7(0;) + 0(x)}„, a t (x) -> on(x) + 9(x) - 6{x (3) 

The parameters appearing in the lattice action differ from the continuum ones, 
but the relation can he calculated exactly with a 2-loop computation in lattice 
perturbation theory 113. rmn 
To find the vortices, we define for each link the quantity □'EMU! 

Y (x,x+i) = + l(x + i) - 7(a0h ~ oti(x). (4) 

Taking a sum around a closed curve C gives the winding number nc- 

Y c = J2 Y i= 2 ™c- (5) 

zee 
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Figure 1: a) The phase diagram of the system, b) The histogram of the volume averages 
of the winding number of a curve of constant size in physical units at various values of the 
lattice spacing a = 1/ (3g ■ In the continuum limit @q — > oo the distance of the peaks seems 
to approach a finite value, which is expected to be independent of the regulator. 



The winding number is a gauge-invariant integer and gives the number of 
vortices going through the curve C 0. Using the difference of only the phase 
angles in Eq. (JiJ) would lead to a non-invariant quantity. 

At the mean-field level the notion of a vortex makes only sense in the Higgs 
phase, but our definition (|J) is perfectly valid in all phases and agrees with 
the intuitive picture of a vortex whenever the bare tree-level potential has a 
degenerate minimum. We calculate the full path integral numerically using 
lattice Monte Carlo simulations, which means including the effect of thermal 
fluctuations to the mean-field picture. 

The quantity we are mainly interested in at this stage is the density of 
thermally generated vortices. The naive way to calculate it is to take the 
absolute value of the winding number of a single plaquette and measure its 
expectation value. However, it turns out that in the continuum limit, this 
quantity approaches a universal quantity w 0.2, which is independent of the 
parameters of the theoryU. The reason is that the winding number of a single 
plaquette is an ultraviolet quantity, and the ultraviolet behavior of the theory 
is given by a free massless complex scalar field. 

To obtain a physical quantity, we take a square curve C and keep its 

a The standard method of locatirtgnvortices by finding the zeros of the Higgs ficldQ has often 
been used also in gauge theories Ej. However, on a lattice one can get rid of all the zeros by 
choosing e.g. the unitary gauge. Even if some other gauge is used, the physical interpretation 
of the vortices found this way is ambiguous. 
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size constant in physical units as we approach the continuum limit. While this 
quantity contains a lot of ultraviolet noise, there is also a physical contribution, 
and the problem is to extract it. An analogous quantity is (<p*(p): it diverges in 
the continuum limit, but the divergent term is constant and can be calculated 
exactly E9. The remaining finite part can then be used to probe the phase 
diagram of the theory. For (|nc|) we cannot subtract the divergence exactly, 
but if it is analytical in the parameters x and y as we expect, the difference 
in (|nc|) above and below the transition line is a physical quantityEI. Some 
numerical evidence for that is shown in Fig. [j]b. 

In practice, it is rather difficult to extract the physical contents of the 
vortex density defined here. However, it shows that in the continuum limit, 
more care is needed than simply interpreting the plaquettes with non-zero 
winding numbers as physical vortices. One has to come up with new observ- 
ables that give better understanding of the problem before a non-perturbative 
understanding of defect formation in gauge theories can be obtained. 
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